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Abstract 

We evaluate the phase-coherent transport of electrons along linear structures of 
varying length, which are made from two types of potential wells set in either a 
periodic or a Fibonacci quasi-periodic sequence. The array is described by a tight- 
binding Hamiltonian and is reduced to an effective dimer by means of a decimation- 
renormalization method, extended to allow for connection to external metallic leads, 
and the transmission coefficient is evaluated in a T-matrix scattering approach. 
Parallel behaviors are found for the energy dependence of the density of electron 
states and of the transmittivity of the array. In particular, we explicitly show that 
on increasing its length the periodic array undergoes a metal-insulator transition 
near single occupancy per dot, whereas prominent pseudo-gaps emerge away from 
the band center in the Fibonacci-ordered array. 
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1 Introduction 



The study of electron transport in mesoscopic systems is a fundamental prob- 
lem in nanophysics and nanotechnology [1,2] and rests on Landauer's idea 
that conductance is proportional to phase-coherent electron transmittivity 
[3]. Within this context a great deal of attention has been given in recent 
years to electronic transport in a variety of structures including polymeric 
and molecular wires, nanotubes and quantum wires, quantum dots, and ar- 
rays and networks of quantum dots. 

A semiconductor quantum dot (QD) is often described as an artificial atom, 
in which an external potential replaces the attractive potential of the nucleus 
to confine charge carriers in all three spatial directions [4,5]. The dot size is of 
the order of the Fermi wavelength in the host material, typically between 10 
nm and 1/im, and the confinement is usually achieved by electrical gating of a 
two-dimensional electron gas, often in combinations with etching of the mate- 
rial. Precise control of the number of electrons in the conduction band of a QD 
has been achieved in GaAs heterostructures [6]. An all-electrically controlled 
quantum dot array (QDA) can be realized by electrodes confining single elec- 
trons to the dot regions, with the electron tunnel between neighboring dots 
being controlled by electrical gating [7]. 

As in the case of real clusters, electronic transport through QDA's can be 
treated theoretically within a tight-binding framework [8]. The tight-binding 
model has been used to study Coulomb effects on QD transmittance [9] , mag- 
netoconductance in chaotic arrays [10], Kondo resonances and Fano antires- 
onances in transport [11], and transport through two-dimensional networks 
[12,13,14]. A simplifying assumption treats a QD as a potential well, thus 
omitting a detailed account of its internal structure. This allows analytical 
results to be obtained for transport properties of QDA's (see for instance 
[15,16]). 

In this work we use the simple tight-binding approach to derive an analyti- 
cal treatment of phase-coherent transport through a QDA model composed 
of two different sets of QD's, as originally proposed by Mardaani and Esfar- 
jani [15]. We treat both a periodic and quasi-periodic arrangement of the two 
types of QD's in the array, and report comparative numerical illustrations 
of the density-of-states (DOS) and of the transmittivity of these arrays for 
a number N of QD's up to 100. Our analytical results allow us to evaluate 
very long arrays without any special numerical effort. Specifically, we use a 
decimation-renormalization method [17,18] to reduce the QDA to an effective 
dimer with renormalized site energies and hopping interactions, and supple- 
ment the method with a generalized Kirkman-Pendry relation [19,20] to ac- 
count for the connection of the QDA to external metallic leads. The T-matrix 
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Fig. 1. Schematic Hamiltonian for (a) a periodic QDA and (b) a Fibonacci QDA, 
for the case N = 12. 

formalism introduced in this paper for the calculation of the QDA transmit- 
tivity is equivalent to the formalism given in Ref. [21] for out-of-equilibrium 
leads in the special case of an infinitesimal bias and is suitable to deal with 
renormalized effective Hamiltonians describing very long chains. 



2 Hamiltonian and density-of-states 



We use a one-dimensional (ID) tight-binding Hamiltonian to describe a QDA 
of equally spaced potential wells. The Hamiltonian can be written as 

N 

H QD = E I l ) I + (*M+i I + !| + * + I )} > (!) 
i=i 



where the site energy e« corresponds to the energy level for an electron in the 
i-th well and = tj+i^ is the hopping energy between the i-th and the 

(i + l)-th well. The connection to external leds will be added later. 

We focus on two types of QDA's. The first is a periodic array in which two 
different types of QD's, A and B say, alternate. In the second the two kinds of 
QD's are arranged according to the Fibonacci sequence, i.e., if we define the 
first QD as Fi = A and the second as F 2 = B, the rest of the chain is built 
by the rule F n = F n _i F„_ 2 . The schematic Hamiltonians for these QDA's are 
shown in Fig. 1. 

The related Green's function at energy E is 

G{E) = 1 . (2) 



Here and below, we consider a complex energy variable defined as E + i r] in 
the limit of vanishing positive imaginary part. To evaluate the DOS D{E) we 
use the Kirkman-Pendry relation [19], which relates the DOS to the matrix 
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element G^n{E) expressing coherence between the first and the last site of 
the QDA, 

D(E) = ±-lm-^\nG 1;N (E). (3) 



The matrix element G\^{E) can be evaluated by reducing the QDA to an 
effective dimer through decimation of the intermediate QD's (see for instance 
[17,18]). The renormalized array contains just the first and the last site and 
its Hamiltonian is expressed as a 2 x 2 matrix, 



H, 



f e[ N - 2 \E) hAE) " 
K t NA (E) e^- 2 \E)j 



(4) 



Here, the effective site and hopping energies can be obtained by the recursive 
relations 



e^(E)=et 1 \E)+t w (E) * t j+hj+2 

E-ef +1 >(E) 

£j%{E) = e j+2 + tj+2,j+i — — r^— iy— tj+i,i(E) , 



t 1 j +2 (E)=t 1J+1 (E) 



E-e^'(E) 
1 



E-e^'(E) 



(5) 
(6) 
(7) 



and tj + i t i = iij+i for j > 1, the initial values being given by the Hamiltonian 
parameters namely Ei°\e) = ei and ti^(E) = t\,2- By direct inversion of 
(E — HQjy(E)) we obtain 



Gi t N(E) 



t N ,i(E) 



[E - e[ N - 2 \E)][E - ?r Z, (E)\ - [hAE)} 2 



(JV-2), 



(8) 



We then use Eq. (3) to calculate the DOS. 

In the illustrative calculations that we report here and in the following sections 
we have used the values of the parameters = —0.25 eV, = 0.25 eV, 
tA,B = 1-1 eV and t B)B = 1.0 eV. In Fig. 2 and Fig. 3 we show the DOS 
for QDA's made from various numbers of QD's in a periodic and a Fibonacci 
configuration. For small values of iV (see Fig. 2) the DOS shows iV peaks 
distributed over the energy range. For long chains with N = 1000 (see Fig. 3) 
the DOS is similar to that of an infinite array: in the periodic sequence (left 
panel in Fig. 3) there are two sub-bands separated by an energy gap of width 
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Fig. 2. DOS of periodic array (left) and Fibonacci-ordered array (right) for the case 
N = 4 (solid line) and N = 20 (dotted line), as a function of energy E referred to 
the center of the spectrum at energy Eq and for a spectral width of 4to. 
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Fig. 3. DOS as in Fig. 2, for the case N = 1000. 

\&a — while in the case of a Fibonacci array (right panel in Fig. 3) the 
quasi-periodicity induces a fragmentation of the spectrum and the appearance 
of pseudo-gaps. 

2. 1 The effects of the leads 



The QDA is next connected to an incoming (I) and an outgoing (r) metallic 
lead, which are described by two additional terms in the Hamiltonian, 



H L,1= E {E \n)(n\+t (\n-l)(n\+c.c.)} 



(9) 



and 



H Lr= J2 {E \n)(n\+t (\n)(n + l\+c.c.)}. 



(10) 



n=N 
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Here the site energy Eq in the leads and the hopping energy to have been 
chosen equal to the center of the spectrum of the QDA and to one fourth of 
its spectral width. The presence of the leads modifies the DOS of the system 
and thus affects the electronic transport through the QDA. We accordingly 
have to decimate the leads and to consistently renormalize the energies of sites 
1 and N in the QDA. 

Equation (3) does not apply in general and in particular, when sites 1 and N 
in the array are not edge sites, a modified Kirkman-Pendry relation as derived 
by Farchioni et al. [20] must instead be used. This relation is 



D(E) = ±Im ^-lnG ljN (E,X) 

71 OA 



(11) 



A=0 



where the real parameter A is introduced to select the sites on which the 
electronic states are counted. In our specific case G\ ]\[(E, A) can be written as 



Gi,n{E, A) 



t N ,i(E) 



[E - e[ N - 2 \E) - £(E)][E - e^~' Z) (E) - £{E)\ - [t hN (E)} 



;i2) 



where E = E + A and the term £{E) = \{E - E ) - [\{E - E ) 2 - tg] 1 / 2 is 
the renormalized contribution of the leads to the energy of sites 1 and N. 

The DOS of a periodic QDA connected to an incoming and an outgoing lead is 
shown in Fig. 4. For small numbers of QD's (N = 4 and N = 12) the effect of 
the leads is dominant and the DOS resembles that of a ID monatomic lattice. 
A depression appears at the center of the spectrum for N = 20, and for a long 
array (N = 100) the effect of the leads becomes irrelevant and the DOS is very 
similar to that shown in the left panel of Fig. 3 for the isolated QDA. Similar 
effects on the DOS of the Fibonacci-ordered QDA are seen in Fig. 5: at low iV 
the peak structure shown in Fig. 2 is washed out by the presence of the leads, 
while at large N the typical features due to quasi-periodicity appear. 



At half filling the Fermi level is located near (E — E )/t = for all config- 
urations, as is indicated in Figs. 4 and 5 by an arrow. In this case the DOS 
at the Fermi level vanishes on increasing the length of the periodic QDA. In 
a Fibonacci-ordered array, instead, the number of states at the Fermi level 
remains finite. 
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Fig. 4. DOS of periodic QDA connected to an incoming and an outgoing lead, as a 
function of energy E referred to the band center at energy Eq and for a bandwidth 
of 4t - Left panel: N = 4 (solid line), TV = 12 (solid line and dots) and N = 20 
(dashed line). Right panel: N = 100. 
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Fig. 5. DOS as in Fig. 4, for a Fibonacci-ordered QDA connected to an incoming 
and an outgoing lead. 



3 Transmission through the array 



We rewrite the Hamiltonian H = -f^Qpj + #l 1 + r °f ^ ne effective dimer 
representing the QDA connected to the leads as the sum of two terms, H = 
Ho + Hi. The first term Hq describes an infinite perfect chain with spacing a, 
site energy Eq and hopping energy to, whose elements (2, 3, . . . , N — 1) have 
been decimated as previously described, and is given by 

H = H L l + H L l + (E - Eo) (|1)(1| + {\N)(N\) + t (\l}(N\ + \N)(1\) .(13) 



The remainder 

H \ = ^QD " {£o(|l)(l| + (\N)(N\) + t (\l)(N\ + \N)(1\)} (14) 
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will be viewed as a "perturbation" determining scattering of incoming waves. 
In these equations the quantities E and to are obtained from Eqs. (5) and (7) 
for j = N — 2 by taking = E and t iji+ i = t for all i. 

The wavefunction \ f) at energy E in the continuous spectrum of H is obtained 
from the wavefunction \k) of the unperturbed periodic Hamiltonian Hq, the 
unperturbed Green's function G°(E) = (E — Hq)^ 1 , and the T matrix T(E) = 
H^l-G ^)' 1 as 

\ v } = \k} + G°T\k), (15) 



T as well as -f/j being 2x2 matrices in the space spanned by |0) and \N). The 
projection of \ip) onto the localized function \n) is 

(n\<p) = (n\k) + £ Gn,iTi, m (m\k) , (16) 

/,m=0,l 



where G° nl = (n\G°\l), T hm = (l\T\m), and (m\k) = e ikma . The expressions 
for these matrix elements are given in [18]. We can then write Eq. (16) in the 
form 

(n\<p) = e thna + (G° NA T ljN + Gl N T NA e- 2ik ^ a + G° NtN T N , N + G^T^) e lkna (17) 
and define the transmittance r and the reflectance p by writing 

(n\<p) = < 



Te ikna (n > 1) 

e ikna + pe -ifena ^ < 



We thus obtain 

t = 1 + G° ni Ti^n + G° hN T NA e~ 2lk{N ~ 1)a + G° NN T NjN + G° hl T 1:1 . (18) 

Notice the difference between this Eq. (18), which is valid for both real and 
effective dimers, and Eq. (8) in Ref. [18] applies only to real dimers. 

The transmission coefficient T is given by T = |r| 2 . Numerical results for 
the periodic and the Fibonacci-ordered QDA are shown in Figs. 6 and 7. The 
features of the transmittivity as a function of the energy parallel those of the 
DOS in Figs. 4 and 5, where the effects of the leads have been included: the 
transmittivity drops as the DOS decreases and vanishes in the energy gaps. 
At half filling in the periodic array the transmittivity at the Fermi level is 
high for very short arrays but rapidly drops on increasing the number of dots, 
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Fig. 6. Transmittivity coefficient for a periodic QDA as a function of the energy E, 
referred to the band center at energy Eq and for a bandwidth of 4to- Left panel: 
N = 4 (solid line), N = 12 (solid line and dots) and N = 20 (dashed line). Right 
panel: N = 100. 
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Fig. 7. Transmittivity coefficient as in Fig. 6, for a Fibonacci-ordered array. 

so that the array becomes insulating. On the contrary, in a Fibonacci-ordered 
array the transmittivity near the Fermi level remains high, so that metallic- 
like conduction is preserved independently of the length of the array. There 
are, however, strong minima in the transmittivity in correspondence to the 
pseudo-gaps in the DOS. 



4 Concluding remarks 

In summary, we have presented a fully analytical method for the calculation 
of the density of electron states and the phase- coherent electron transport co- 
efficient in a linear array made of an arbitrary number of two different types 
of potential well with different energy levels. We have specifically considered 
both a periodic sequence of pairs of potential wells and a sequence in which 
the two types of potential wells are arranged according to the Fibonacci se- 
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quence. In both cases the connection of the array to external metallic wires has 
been included. The method does not explicitly diagonalize the tight-binding 
Hamiltonian of the system, but uses a decimation-renormalization procedure 
to reduce the array to an effective dimer, which is then treated by a suitably 
adapted T-matrix scattering approach. 

The density of states and the transmittivity vanish at the same energy when 
the connection to external metallic wires is included. The calculations ex- 
plicitly show that an energy gap opens in the periodic array as its length is 
increased, so that metallic conductance turns into insulating behavior in the 
case of single occupancy of each well. With increasing length the Fibonacci- 
ordered array develops instead pseudo-gaps away from the band center, and 
again a metal-to-insulator transition is expected for appropriate values of the 
filling factor. 
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